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In this paper we address the question how the control of delayed measured chaotic systems can be 
improved. Both unmodified OGY control and difference control can be successfully applied only for 
a certain range of Lyapunov numbers depending on the delay time. We show that this limitation can 
be overcome by at least two classes of methods, by rhythmic control and by the memory methods 
linear predictive logging control (LPLC) and memory difference control (MDC). 
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I. INTRODUCTION 

Delay is a generic problem in the control of chaotic 
systems. The effective delay time r in any feedback loop 
is the sum of at least three delay times, the duration 
of measurement, the time needed to compute the appro- 
priate control amplitude, and the response time of the 
system to the applied control. The latter effect appears 
especially when the applied control additionally has to 
propagate through the system. These response time may 
extend to one or more cycle lengths Q. 

For the formal situation of fixed point stabilization in 
time-continuous control, the issue of delay has been in- 
vestigated widely in control theory, dating back at least 
to the Smith predictor 0. This approach mimics the, 
yet unknown, actual system state by a linear prediction 
based on the last measurement. Its time-discrete counter- 
part (LPLC) discussed below allows to place all eigenval- 
ues of the associated linear dynamics to zero, and always 
ensures stability. The (time-continuous) Smith predictor 
with its infinite-dimensional initial condition had to be 
refined 0, , giving rise to the recently active fields of 
model predictive control |jj . For fixed point stabilization, 
an extension of permissible latency has been found for a 
modified proportional-plus-derivative controller HJ. 

If one wants to stabilize the dynamics of a chaotic sys- 
tem onto an unstable periodic orbit, one is in a special sit- 
uation. In principle, a proper engineering approach could 
be to use the concept of sliding mode control Jj|, i. e. to 
use a co-moving coordinate system and perform suitable 
control methods within it. However, this requires the 
quite accurate knowledge of the whole trajectory (and 
the direction of the stable manifold) with the respective 
numerical or experimental costs. 

Therefore direct approaches have been developed by 
cxplicitely taking into account either a Poincare sur- 
face of section [8( or the explicit periodic orbit length 
0- This field of controlling chaos, or stabilization of 
chaotic systems, by small perturbations, in system vari- 
ables |ldj| or control parameters Q , emerged to a widely 
discussed topic with applications in a broad area from 
technical to biological systems. Especially in fast systems 
[Til or for slow drift in parameters 0, 0], differ- 
ence control methods have been successful, namely the 
time-continuous Pyragas sche me pjl , ETDAS an d 
time-discrete difference control [l5| . 



Like for the control method itself, the discussion of 
the measurement delay problem in chaos control has to 
take into account the special issues of the situation: In 
classical control applications one always tries to keep the 
control loop latency as short as possible. In chaotic sys- 
tems however, one wants to control a fixed point of the 
Poincare iteration and thus has to wait until the next 
crossing of the Poincare surface of section, where the sys- 
tem again is in vicinity of that fixed point. 

The stability theory and the delay influence for time- 
continuous chaos control schemes has been studied exten- 
sively 0, 0, 0, 0, |2(j , and an improvement of control 
by periodic modulation has been proposed in |2l|. For 
measurement delays that extend to a full period, how- 
ever no extension of the time-continuous Pyragas scheme 
is available. 

In this paper we investigate time-discrete control 
schemes and focus on the question what limitations occur 
if one applies OGY control @ or difference feedback 0] 
in the presence of time delay, and what strategies can 
be used to overcome these limitations. We show how the 
measurement delay problem can be solved systematically 
for OGY control and difference control by rhythmic con- 
trol and a memory method and give constructive direct 
and elegant formulas for the deadbeat control in the time- 
discrete Poincare iteration. While the predictive control 
method LPLC presented below for OGY control has a 
direct correspondence to the Smith predictor and thus 
can be reviewed as its somehow straightforward imple- 
mentation within the unstable subspace of the Poincare 
iteration, this prediction approach does not guarantee a 
stable controller for difference control. However, within a 
class of feedback schemes linear in system parameters and 
system variable, there is always a unique scheme where 
all eigenvalues are zero, the MDC scheme presented be- 
low. The method can be applied also for more than one 
positive Ljapunov exponent, and shows, within validity 
of the linearization in vicinity of the orbit, to be free 
of principal limitations in Ljapunov exponents or delay 
time. For zero delay (but the inherent 1 period delay of 
MDC), MDC has been demonstrated experimentally for 
a chaotic electronic circuit [T^ | and a thermionic plasma 
discharge diode 0], with excellent agreement, both of 
stability areas and transient Ljapunov exponents, to the 
theory presented here. 

The paper is organized as follows. After introducing 
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the notation within a recall of OGY control, we give 
a brief summary what limitations occur for unmodified 
OGY control; details can be found in j^. In Section ITTT1 
we introduce different memory methods to improve con- 
trol, of which the LPLC approach appears to be superior 
as it allows stabilization of arbitrary fixed points for any 
given delay. The stabilization of unknown fixed points 
is discussed in Section IIVI where we present a memory 
method (MDC) that again allows stabilization of arbi- 
trary unstable fixed points. As shown in Appendix 
for all systems with only one instable Lyapunov number, 
the iterated dynamics can be transformed on an eigen- 
system which reduces to the one-dimensional case. The 
explicit formulas for the case of higher-dimensional sub- 
spaces are given in Appendices |5] and |0 



II. CONTROL OF UNSTABLE PERIODIC 
ORBITS 

A. Ott, Grebogi, Yorke control 

The OGY method given by Ott, Grebogi and Yorke 
8] stabilizes unstable fixed points (or unstable periodic 
orbits utilizing a Poincare surface of section) by feedback 
that is applied in vicinity of the fixed point x* of a dis- 
crete dynamics x t +i = f{x t ,r). 

Thus, for a chaotic flow (or corresponding experiment) 
one reduces the system dynamics 



F(x,r) 



(1) 



to the discrete dynamics between subsequent Poincare 
sections at to,ti,...t n . This description is fundamen- 
tally different from a stroboscopic sampling as long as 
the system is not on a periodic orbit, where the sequence 
of differences (ti — ii_i) would show a periodic structure. 

If there is only one positive Ljapunov exponent, we 
can proceed considering the motion in unstable direction 
only (see Appendix 0, i. e. a one-dimensional iterated 
map. (For 2 or more positive Ljapunov exponents on can 
proceed in a similar fashion, see Add. IBlandlGll 

In OGY control, the control parameter r t is made time- 
dependent. The amplitude of the feedback r t = r — tq 
added to the control parameter ro is proportional by a 
constant e to the distance x — x* from the fixed point, 
i. e. r = ro + e(xt — x*), and the feedback gain can be 
determined from a linearization around the fixed point, 
which reads, if we neglect higher order terms, 

fdf 

f{x t ,r + r t ) = f(x*,r ) + (x t - x) ■ I — 

• (§0 , 

= f(x*,r Q ) + \(x t ~x*) + nr t 

= f(x*,r ) + (\ + ne)-(x t -x*) (2) 

The second expression vanishes for e = that is, 

in linear approximation the system arrives at the fixed 



point at the next time step, x t +\ = x*. The uncontrolled 
system is assumed to be unstable in the fixed point, i. e. 
|A| > I. The system with applied control is stable if the 
absolute value of the eigenvalues of the iterated map is 
smaller than one, 



Ft+i 



'I = |(A + fie) ■ (x t - x*)\ < \x t 



(3) 



Therefore e has to be chosen between (—1 — X)//j, and 
(+1 — A)//x, and this interval is of width 1j\x and inde- 
pendent of A, i.e. fixed points with arbitrary A can be 
stabilized. This property however does not survive for 
delayed measurement. One has to develop further con- 
trol strategies to cover this class of systems. 



B. Delay matching in experimental situations 

Before discussing the time-discrete reduced dynamics 
in the Poincare iteration, it should be clarified how this 
relates to an experimental control situation. On the first 
glance, the time-discrete viewpoint seems to correspond 
only to a case where the delay (plus waiting time to the 
next Poincare section) exactly matches the orbit length, 
or a multiple of it. The generic experimental situation 
however comes up with a non-matching delay. Appli- 
cation of all control methods discussed here requires to 
introduce an additional delay, usually by waiting for the 
next Poincare crossing, so that measurement and control 
are applied without phase shift at the same position of 
the orbit. In this case the next Poincare crossing po- 
sition Xt+i is a function of the values of x and r at a 
finite number of previous Poincare crossings only, i. e. 
it does not depend on intermediate positions. Therefore 
the (a priori infinite-dimensional) delay system reduces 
to a finite-dimensional iterated map. 

If the delay (plus the time of the waiting mechanism 
to the next Poincare crossing) is not matching the orbit 
length, the control schemes may perform less efficient. 
Even for larger deviations from the orbit, the time be- 
tween the Poincare crossings will vary only marginally, 
thus a control amplitude should be available in time. In 
practical situations therefore the delay should not exceed 
the orbit length minus the variance of the orbit length 
that appears in the respective system and control setup. 

In a formal sense, the Poincare approach ensures ro- 
bustness with respect to uncertainties in the orbit length, 
as it always ensures a synchronized reset of both trajec- 
tories and control. Between the Poincare crossings the 
control parameter is constant, the system is independent 
of everything in advance of the last Poincare crossing. It 
is solely determined by the differential equation (or ex- 
perimental dynamics). Thus the next crossing position 
is a well defined iterated function of the previous one. 

This is quite in contrast to the situation of a delay- 
differential equation (as in Pyragas control) , which has an 
infinite-dimensional initial condition it 'never gets rid of. 
One may proceed to stability analysis via Floquet theory 
[25| as investgated for continuous and Poincare-based 
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23, 26] control schemes. Though a Poincare crossing de- 
tection may be applied as well, the position will depend 
not only on the last crossing, but also on all values of 
the system variable within a time horizon defined by the 
maximum of the delay length and the (maximal) time 
difference between two Poincare crossings (due to the 
non-stroboscopic character). Thus the Poincare itera- 
tion would be a function between two infinite-dynamical 
spaces. Apart from further mathematical subtleties, for a 
delay differential equation with fixed delay the major ad- 
vantage of a Poincare map, reducing the system dynam- 
ics to a low-dimensional system, thus completely breaks 
down. 

For all control schemes discussed within this paper, 
howvever the additional dimensionality is not a contin- 
uous horizon of states, but merely a finite set of values 
that were measured at the previous Poincare crossings. 



III. OGY CONTROL: DELAYED CASE, 
RHYTHMIC AND SIMPLE MEMORY CONTROL 

A. Unmodified control of delayed measured maps 

To illustrate the problem, we recall the simple case 
where no modification of the OGY scheme is taken into 
account [22|. For r time steps delay, unmodified propor- 
tional feedback is applied: 



r t = e(x t -r - x*). 



(4) 



Without loss of generality, we can choose rt = if no 
control is applied, and x* = in the remainder. Using 
the time-delayed coordinates (xt, Xt-i, x t ~2, ■ ■ ■ xt- T ) T '- 



* f 



eA*\ / xt \ 
' x 



1 



\x t -T+i/ \0 ■ 



(5) 



• 010/ \xt-J 



To ensure stability, one has to fulfill \oii\ < 1 for all eigen- 
values ccj . The stability area is investigated in [22| , yield- 
ing a r-dependent maximal Ljapunov number of 



An 



f 

f + -. 

r 



(6) 



For A > A max (r) no control is possible |22||. 



B. Rhythmic OGY control 

As pointed out for difference control in the case r = 
in [13, 123, |2Jj u3 , one can eliminate the additional de- 
grees of freedom caused by the delay term. One can 



restrict himself to apply control rhythmically only ev- 
ery t + 1 timestcps (t + 2 for difference control), and 
then leave the system uncontrolled for the remaining 
timesteps. Then e = e(t) appears to be time-dependent 
with 



e(t mod t) = (e , 0, . . . , 0) 



(7) 



and, after (t+1) iterations of JSJ, we again have a matrix 
as in J5J), but with A T+1 instead of A. Equivalently, we 
can write 



x t+ {r+i) = A T+1 x t + e pLX t . 



(8) 



What we have done here, is: controlling the (r + 1)- 
fold iterate of the original system. This appears to be 
formally elegant, but leads to practically uncontrollable 
high effective Lyapunov numbers A r+1 for both large A 
and large r. 

Even if the rhythmic control method is of striking sim- 
plicity, it remains unsatisfying that control is kept quiet, 
or inactive, for r time steps. Even if the state of the 
system x is known delayed by r, one knows (in principle) 
the values of xt for t < r, and one could (in principle) 
store the values 5r t ~ T ■ ■ ■ 5r t of the control amplitudes 
applied to the system. This can be done, depending on 
the timescale, by analog or digital delay lines, or by stor- 
ing the values in a computer or signal processor (there 
are some intermediate frequency ranges where an exper- 
imental setup is difficult). 

Both methods, rhythmic control and simple feedback 
control in every time step, have their disadvantages. 

For rhythmic control it is necessary to use rather large 
control amplitudes, in average A T /r, and noise sums up 
to an amplitude increased by factor y/r. 

For simple feedback control the dimension of the sys- 
tem is increased and the maximal controllable Lyapunov 
number is bounded by ©. 

One might wonder if there are control strategies that 
avoid these limitations. This has necessarily to be done 
by applying control in each time step, but with using 
knowledge what control has been applied between the 
last measured time step t — r and t. This concept can 
be implemented in at least two ways, by storing previous 
values of Xt or previous values of 5r t (LPLC). 



C. Control using memory for previous states 

The first memory method extends the single delay line 
to several delay lines, using one control gain coefficient 
for each: 

r t = eiXt-i + e 2 x t -2 + ■ ■ ■ + s n+1 x t -„-i (9) 
For n steps memory (and one step delay) the control 
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matrix is 



/Xt+l\ A El ••• 

' 1 1 

i '•. 



£n En+l\ / Xt \ 

4 



(10) 



\Xt-J \0 1 / \x t -n-J 

with the characteristic polynomial 

n 

(a- X)a n+1 +^ £l a™- 



(11) 



We can choose eci = = .. . a„+2 = — A/(n + 2) and 
evaluate optimal values for all £j by comparing with the 
coefficients of the product YIa=i i a ~ a i)- This method 
allows control up to A max = 2 + n, therefore arbitrary A 
can be controlled if a memory length of n > A — 2 and 
the optimal coefficents e% are used. 

For more than one step delay, one has the situation 
El = 0, . . . , £ T -x = 0. This prohibits the 'trivial pole 
placement' given above, (choosing all on to the same 
value) and therefore reduces the maximal controllable A 
and no general scheme for optimal selection of the £j 
applies. One can alternatively use the LPLC method de- 
scribed below, which provides an optimal control scheme. 
One could wonder why to consider the previous state 
memory scheme at all when it does not allow to make all 
eigenvalues zero in any case. First, the case of up to one 
orbit delay and moderately small A already covers many 
low-period orbits. Second, there may be experimental 
setups where the feedback of previous states through ad- 
ditional delay elements and an analog circuit is experi- 
mentally more feasible than feedback of controlled am- 
plitudes. 

In summary, all three methods discussed in this sec- 
tion have principal restrictions, but may be of advantage 
in special situations especially when a simple setup is 
required. 



D. Linear predictive logging control (LPLC) 

If it is possible to store the previously applied control 
amplitudes r t ,r t _i, . . ., then one can predict the actual 
state Xt of the system using the linear approximation 
around the fixed point. That is, from the last measured 
value Xt-r and the control amplitudes we compute esti- 
mated values iteratively by 



Vt-i+i = Aa; t _i + fin-t 



(12) 



leading to a predicted value yt of the actual system state. 
Then the original OGY formula can be applied, i. e. r t — 
—\/fiy t . In this method the gain parameters are again 



linear in Xt- T and all {r t >} with t — r < t' < t, and the 
optimal gain parameters can be expressed in terms of A 
and \x. 

In contrast to the memory method presented in the 
previous subsection, the LPLC method directs the sys- 
tem (in linear approximation) in one time step onto the 
fixed point. However, when this control algorithm is 
switched on, one has no control applied between t — t 
and t — 1, so the trajectory has to be fairly near to the 
orbit (in an interval with a length of order 5/X r , where 
6 is the interval halfwidth where control is switched on). 
Therefore the time one has to wait until the control can 
be successfully activated is of order A T_1 larger than in 
the case of undelayed control. 

The LPLC method can also be derived as a general 
linear feedback in the last measured system state and 
all applied control amplitudes since the system was mea- 
sured; choosing the feedback gain parameters so that the 
linearized system has all eigenvalues zero. The linear 
ansatz 



e ■ xt-r-i + rji r t-i + ■■■ Vrn-i 



(13) 



leads to the dynamics in combined delayed coordinates 

(x t ,x t - 1, ... ,x t - T ,rt-i, ■ ■ . ,n- T ) 



I xt+1 \ /A 



Xt-T + l 



e jji Tf2 







1 

e rj\ r\2 
1 



'It 



\r t - T+ J V 

giving the characteristic polynomial 



Xt-l 



1 



= -a T {a T+1 +a T {-\-r ]1 ) 

+ a T ^ 1 (X ■ 771 - r] 2 ) + a T ~ 2 (\ ■ 7/2 - m) 
. . . + a : (A • T7 T _i - ri T ) + (A • Vt - e)). (14) 

All eigenvalues can be set to zero using s = — A r+1 and 
rji = — X 1 . The general formulas even for more than one 
positive Lyapunov exponent or multiparameter control 
are given in Appendix iBl 



E. Nonlinear predictive logging control 

One can also consider a nonlinear predictive logging 
control (NLPLC) strategy as the straightforward exten- 
sion to the LPLC method for nonlinear prediction. If 
the system has a delay of several time steps, the interval 
where control is achieved becomes too small. However, if 
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it is possible to extract the first nonlinearities from the 
time series, prediction (and control) can be fundamen- 
tally improved. In NLPLC the behaviour of the system 
is predicted each time step by 



x t +i 



A 2 



f*2 2 / 1 \ 

+ M r * + -7r r t + ux tn + o(x t , . . .) 



with using the applied control amplitudes {r t } for each 

time step. This equation has to be solved for rt using 
i 

x t +i = 0. 

A similar nonlinear prediction method has been de- 
scribed by Petrov and Showalter [2j| . They approximate 
the Xt+i(xt,r t ) surface directly from the time series and 
use it to direct the system to any desired point. Both 
Taylor approximation or Petrov and Showalter method 
can be used here iteratively, provided one knows the de- 
lay length. 

Both approaches could be regarded as a nonlinear 
method of model predictive control 5], applied to the 
Poincare iteration dynamics. 

From a practical point of view, it has to be mentioned 
that one has to know the fixed point x* more accurate 
than in the linear case. Otherwise one experiences a 
smaller range of stability and additionally a permanent 
nonvanishing control amplitude will remain. This may be 
of disadvantage especially if the fixed point drifts in time 
(e.g. by other external parameters such as temperature) 
or if the time series used to determine the parameters is 
too short. 



is given by 



' ] 1 

1 



en -e/A ( Xt \ 
' x 



\xt-J \0 ■■■ 











(15) 



J \Xt-T-J 



in delayed coordinates yj = {xt,xt~i, ■ ■ ■ x t - T -\)- As we 
have to use Xt- T -i in addition to x t - Tl the system is 
of dimension r + 2, and the lower bound of Lyapunov 
numbers that can be controlled arc found to be 1221 



3 + 2r 
" 1 + 2t ' 



(16) 



The controllable range is smaller than for unmodified 
OGY control, and is restricted to oscillatory repulsive 
fixed points with Ai„f < A < — 1. A striking observation 
is that inserting T+k for r in eq. © leads exactly to the 
expression in eq. (|16|) which reflects the fact that the dif- 
ference feedback can be interpreted as a discretized first 
time derivative, taken at time t — t — i. For details and 
stability area see p^|. 



B. Rhythmic difference control 



are obtained only for —3 < A < +1, so 



IV. STABILIZATION OF UNKNOWN FIXED 
POINTS 

As all methods mentioned above require the knowl- 
edge of the position of the fixed point, one may wish 
to stabilize purely by feeding back differences of the sys- 
tem variable at different times. Without delay, difference 
feedback can be used successfully for e/i = —A/3, and 
eigenvalues of modulus smaller than unity of the matrix 
A + s/j, —efi 
1 

this method stabilizes only for oscillatory repulsive fixed 
points with -3 < A < -1 [Tll^ . 

Due to the inherent one period delay of MDC, the r 
period delay case of MDC corresponds, in terms of the 
number of degrees of freedom, to the r + 1 period delay 
case of LPLC. 



A. Unmodified difference control of iterated 
delayed measured maps 

In the presence of r steps delay the linearized dynamics 
of a simple difference feedback rt — rg = e(x t - T — x t - T -i) 



To enlarge the range of controllable A, one again has 
the possibility to reduce the dimension of the control pro- 
cess in linear approximation to one by applying control 
every r + 2 time steps. 



x t+i — Axt + iie(xt- T — Xt-r 
= (A T+1 + iie\ — fis)x t - 

and the goal Xt+i = can be fulfilled by 



l) 



(17) 



A T 



[IE 



1 - A 



l-A 



(18) 



to achieve 



One has to choose fie between /i£± = — - 
control as shown in Fig. [T] The case r — has already 
been discussed in 0, lU El E3 • With rhythmic control, 
there is no range limit for A, and even fixed points with 
positive A can be stabilized by this method. 

When using differences for periodic feedback, one still 
has the problem that the control gain increases by A r , 
and noise sums up for t + 1 time steps before the next 
control signal is applied. Additionally, now there is a 
singularity for A = +1 in the "optimal" control gain given 
by (I17|) . This concerns fixed points where differences Xt — 
x t -i when escaping from the fixed point are naturally 
small due to a A near to +1. 
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Here one has to decide between using a large control 
gain (but magnifying noise and finite precision effects) 
or using a small control gain of order /xe_(A = +1) = 
t + i (but having larger eigenvalues and therefore slow 
convergence) . 




FIG. 1: Periodic difference feedback for r = 0,1,2,3,4,5: 
Maximal, optimal and minimal value of jj,e for given A to 
obtain stabilization by control applied every r + 2 time steps. 



Two other strate gies that have been discussed by So- 
colar and Gauthier [29j are discretized versions of time- 
continuous methods. Control between A = — (3 + i?)/(l — 
R) and A = -1 is possible with discrete-ETDAS (R < 1) 
r t = ^Efclo ^(^f-it — Xt-k-i) and control between A = 
— {N + 1) and A = -1 is acheived with discrete- NTDAS 

(let N be a positive integer) r t — e i^Xt — 12k=o x t-k) ■ 
Both methods can be considered to be of advantage even 
in time-discrete control in the Poincare section, e.g. if 
the number of adjustable parameters has to be kept 
small. Whereas these methods are mainly applied in 
time-continuous control, especially in analogue or optical 
experiments, for time-discrete control the MDC strategy 
described below allows to overcome the limitations in the 
Lyapunov number. 



C. Memory difference control (MDC) 

One may wish to generalize the linear predictive feed- 
back to difference feedback. In contrary to the LPLC 
case, the reconstruction of the state Xt- T from differences 
and applied control amplitudes r t _j is 
no longer unique. As a consequence, there are infinitely 
many ways to compute an estimate for the present state 



of the system, but only a subset of these leads to a con- 
troller design ensuring convergence to the fixed point. 

Under these there exists an optimal every-step control 
for difference feedback with minimal eigenvalues and in 
this sense optimal stability. This memory difference con- 
trol (MDC) method has been demonstrated in an elec- 
tronic experiment and a plasma diode [l4| . 

To derive the feedback rule, we directly make the linear 
ansatz 

r t = e • {xt-r-i - xt-r-i-i) + mn-i + ■■■ rh-n-T (19) 
leading to the dynamics in combined delayed coordinates 

(x t + l,x t ,..., x t -r+i,rt-i, ■ • ■ , r t _ T+ i) T 

= M • (x t ,x t - 1, . . .,x t -r,r t -i, . ■ - ,r t -r) T (20) 

with 

/A ■■■ e — e rji rfi - r)i 7) T — J7t+i\ 

1 

1 '■• 



■ ■ ■ £ — £ 7)1 7)2 - r/i Tfr — r) T + i 

1 



V 1 J 

giving the characteristic polynomial 

= -a T (a T+1 +a T (-\- m ) 

+a T_1 (A • r/i - r? 2 ) + a T ~ 2 (A • r\ 2 - r\ z ) 
... + a 2 (\ ■ r\ T - 2 - Vr-l) 
+a 1 (\-T lT - 1 -r ]T -e) + (\-r lT +e)). (21) 

All eigenvalues can be set to zero using e = — A r+1 /(Z — 1) 
r h = +\ T /(l - 1) and m = -X i for 1 < i < t - 1. 
The general formulas even for more than one positive 
Lyapunov exponent or multiparameter control arc given 
in Appendix [HI 

V. CONCLUSIONS 

We have presented methods to improve Poincare- 
section based chaos control for delayed measurement. 
For both classes of algorithms, OGY control and differ- 
ence control, delay affects control, and improved control 
strategies have to be applied. Improved strategies con- 
tain one of the following principle ideas: Rhythmic con- 
trol, control with memory for previous states, or control 
with memory for previously applied control amplitudes. 

In special cases the unmodified control, previous state 
memory control, or rhythmic control methods could be 
considered, especially when experimental conditions re- 
strict the possibilities of designing the control strategy. 
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In general, the proposed LPLC and MDC strategies 
allow a so-called deadbeat control with all eigenvalues 
zero; and they are in this sense optimal control methods. 
All parameters needed for controller design can be cal- 
culated from linearization parameters that can be fitted 
directly from experimental data. This approach has also 
been sucessfully applied in an electronic |l3j and plasma 
|14| experiment. 

APPENDIX A: TRANSFORMATION ON THE 
EIGENSYSTEM 

Here we derive how for one unstable dimension the sta- 
bilization problem reduces to the one-dimensional case. 
Using a covariant basis from right eigenvectors e u and 
e s to eigenvectors A u und A s and the corresponding con- 
travariant left eigenvectors / u and / s of matrix L. we 
can transform the linearized dynamics 

$ t+1 =Lx t + Mr t (Al) 

with help of 1 = e u ■ -f u + e s ■ -f s and L = A u e u • -f u + 
X s e s ■ -f s . We define as coordinates in the eigensystem 
xf := f u x t and x s t := f s x t giving 

%t+i = ^uXt + npu 

x s t+1 = \ s xl + r t fj> s (A2) 

We consider a general ansatz that the control signal is 
linear in x t ~j and r t -i, 

oo oo 

n = E K J ' + E ■ T t-t (A3) 

3=0 i=0 

Here rt = K ■ xt is OGY control. For rt follows 

oo oo 

n = E Kjfa ' 'f u + ^ s ' 'f s )*t-j + E ^ ' n-i 

3=0 i=0 
oo oo oo 

= E^ • 2*)xt-j + E^ ■ e s )x s t _ J +J2 7 lf r *-« 

3=0 j=0 i=0 

oo oo oo 

= E K >u + E K PU + E * ■ r *-i ( A4 ) 

j=o ]=a »=o 

If control in stable direction is chosen zero, i.e. Vjifj = 
0, the dynamics decouples in two two systems in stable 
direction (that is unaffected, but was stable before) and 
in unstable direction, 

oo oo 

x? +1 = \ u x? + Mu (X Kfx^ + E ^ • r ^)- ( A5 ) 

3=0 i=0 

This equation is only one-dimensional, so it is sufficient 
to investigate the one-dimensional control problem if 
there is only one instable direction. The generalization 
to higher-dimensional unstable subspaces is straightfor- 
ward. 



APPENDIX B: STABILIZATION OF DELAYED 
MAPS 

We consider the motion around an unstable fixed point 
x* of the map 

x t+ i =F{x u r t ) (Bl) 

where f t and Xt have the same dimension of the phase 
space of the system (although it is desirable to achieve 
control with a minor number of control parameters). The 
linear time evoution around the instable fixed point (we 
choose x* — and f — in the fixed point) is given by 
Jacobians D x =: L and D r =: M 

x t+ i = Lx t + Mf t (B2) 

Delayed measurement means that the applied control f t 
can only be computed from an older value Xt—r- Addi- 
tionally one can store the control amplitudes f\. . .f%—r 
applied in the meantime to the system and use the lin- 
earization to predict the prese nt p osition x (Linear pre- 
dictive logging control, LPLC |23|). 

Thus we have the general feedback ansatz 

T 

r t =Kx t -T + ^N j r t - j . (B3) 

3=1 

Under the condition that detM ^ 0, the feedback can 
be chosen to 

K = -M _1 I7 +1 (B4) 
Vi<i< T N. L = -M- X VM (B5) 

Iterating x t — L T xt- T + SJ=i I^^Mft—j, one has 
x t +i = Lx t + Mr t = L T+1 x t - T - MM" 1 L T+1 a; t _ T 

T T 

+ X L j - l Mr t -j - X MM- l JJM? t -j = 0. 

3 = 1 3=1 

APPENDIX C: DELAYED MEASURED MAPS - 
STABILIZING UNKNOWN FIXED POINTS 

Now we show that even fixed points whose exact posi- 
tion is not given (only an approximative value is needed 
to determine the position of a <5-ball inside which con- 
trol is switched on) can be stabilized by a difference con- 
trol method that is similar to the LPLC method given 
above, even if one has only delayed knowledge of differ- 
ences of the system variables Xt- T — Xt— T —i- Combined 
with stored values of the meantime control amplitudes 
f\. . .f t -r-i we propose the control scheme 

T + l 

ft = K(S t - T - Xt-r-l) + E N i?t-i ( C1 ) 

3=1 
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with the feedback matrices 

K = -M~ l L r+2 {L - 



(C2) 

Vi<. t < T Ni = -M- X VM (C3) 
iV r+ i = M -1 i T+1 (£- Vf x M (C4) 



Here we have to assume that not only M but also (L — 1) 
is invertible. Using the linear approximation one easily 
computes directly that this control leads to x t +i — 0, and 
again this is a so-called deadbeat control scheme where 
all eigenvalues are zero. 

Using the linearized dynamics, we have 

x t -T = Lxt-r-i + Mft-r-i (C5) 

r+l 

x t = L T+1 x t ^ 1 + Y,L ] - 1 Mf t - J (C6) 
i=i 

and the next iteration reads 

T+l 

x t+1 = Lx t + Mf t = i T+2 f 4 _ r _i + ^ jJMPt-j 



—MM~ 1 L T+2 (L - Vj^LSt-r-i 
-MM- l L T+2 (L - Vj^Mft-r-i 
+MM- 1 L T+2 (L - lyt-xt-r-i 

T 

-^MM- X V Mft-j 
+MM~ 1 L T+1 {L - l^Mft-r-i 



(C7) 



By this method the delay can be overcome and even 
fixed points whose position is known inaccurate can 
be stabilized from measured differences, provided that 
one has a sufficient number of control parameters, i. e. 
dimr = dimx, and both M and [L — 1) are invertible. 
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